Diffusion-limited-aggregation on a directed small world network 
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For real world systems, nonuniform medium is ubiquitous. Therefore, we investigate the diffusion- 
limited-aggregation process on a two dimensional directed small-world network instead of regular 
lattice. The network structure is established by rewiring connections on the two dimensional directed 
lattice. Those rewired edges are controlled by two parameters 6 and m, which characterize the 
spatial length and the density of the long-range connections, respectively. Simulations show that 
there exists a maximum value of the fractal dimension when 6 equals zero. Interestingly, we find 
that the symmetry of the aggregation pattern is broken when rewired connections are long enough, 
which may be an explanation for the formation of asymmetrical fractal in nature. Then, we perform 
multifractal analysis on the patterns further. 

PACS numbers; 



I. INTRODUCTION 



Nonequilibrium growth models leading naturally to 
self-organized fractal patterns, the structure of which 
strongly depends on the dynamics of the growth pro- 
cess, are of continuing interest due to their relevance for 
many important fields 1]. Diffusion-limited -aggregation 
(DLA) 121 is probably the most remarkable growth model 
for pattern formation. This model generates complex 
and mysterious fractal structures[3,l3i which seem to be 
generated as well in natural systems in which growth is 
controlled by diffusive, including dielectric breakdown^], 
electrochemical deposition j^L colloid aggregation^, film 
growth^, viscous fingeringj^, Laplacian flow[Tol| etc. 

In the original DLA modelQ, particles released at a 
point distant from the cluster execute random walks un- 
til they find a nearest neighbor site of the cluster and 
irreversibly stick at this site. However, for many physical 
processes, the medium is nonuniform so that the prob- 
ability of jumping from site i to site j is usually not 
equal to that from j to i. Moreover, except nearest- 
neighboring jumps, there also exist some nonlocal jumps 
through which the particle can move to a distant site 
at a step. One example is the diffusion of adatoms on 
metal surfa ces, in which the long jumps play a dominat- 
ing role [ill A bunch of defects or impurities in the 
substrate may also play the part of the long-jump path 
in the case of weak absorbate-substrate interactionp^. 
which is important for the thin-film growth, heteroge- 
neous catalysis, and oxidation. Hence, the traditional 
undirected regular lattice will miss important informa- 
tion of medium and it is unconformity to characterize 
the actual DLA process. This thus calls for the use of 



network structure with directed long-range connections. 

On the other hand, complex networks have rece ntly 
attracted an increasing interest among phvsicists|l4l llSl 
IT^ . In particular, small-world (SW) networks, intro- 
duced by Watts and Strogatzj^, have been extensively 
studied because they constitute an interesting attempt 
to translate the complex topology of social, biological, 
and physical networks into a simple model. Two dimen- 
sional (2D) SW networks result from randomly rewiring 
a fraction of links of a 2D regular lattice. Several dy- 
namical models have recently been studied in order to 
understand the effect of SW topology on classical systems 
such as the Ising model th e sp read of epidemics [l^. 
aggregation 20J , random walks [2H. etc. Such models are 
expected to capture the essential features of the compli- 
cated processes taking place on real networks. 

In this paper, we investigate the DLA process on a 2D 
directed SW network, in which the directed links corre- 
spond to the directed irreversibly jumps and the node is 
regarded as the lattice point of the real space, respec- 
tively. 



II. THE MODEL 

In order to construct the directed SW network, we 
start from a 2D square lattice of size L x L consisting 
of sites linked to their four nearest neighbors by both 
outgoing and incoming links, as shown in Fig. 1. Then, 
we reconnect nearest-neighbor outgoing links to a differ- 
ent site chosen at random with the probability 
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where r is the lattice distance between the two selected 
sites and 9{0 £ [—1, 1]) is the clustering exponent which 
characterizes the localized degree of the long-range links 
in the lattice. The formula corresponds to the diffusiv- 
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FIG. 1: Sketch map of a directed 2D SW lattice. Dotted lines 
represent rewired connections. Arrows indicate the direction 
of the corresponding connection. 
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FIG. 2: The simulation patterns as a function of the clustering 
exponent for the number of long-range connections m=1000, 
10000 and 100000 respectively. Each color represents 1000 
particles in order. 



ity represented by Arrhenius relation as usual in sur- 
face science [III ITa IT^ . The reconstructing process is 
repeated until m, the number of long-range rewiring con- 
nections, reaches a desired value. Note that by perform- 
ing this procedure every site will have exactly four fixed 
outgoing links and a random number of incoming links. 
When the clustering exponent 9=0, we have a uniform 
distribution over long-range connections, and the present 
model reduces to the basic 2D directed SW network[2^ . 
As 9 ^ 1 {6 ^ 1 denotes 9 tends to 1), the long-range 
links of a site become more and more local in its vicinity. 
In reverse, as — > — 1, the long-range rewiring outgoing 
links are in favor of the farther sites. Thus, the clustering 
exponent 9 serves as a structural parameter controlling 
the spatial length of the long-range connections. 

Based on the directed 2D SW network as constructed 
above, we have performed extensive numerical simula- 
tions for the DLA with size of the reconstructed lat- 
tice L=1000 with number of particles iV=10000. Star- 
ing from an immobile seed at the center of the lattice, a 
particle is released at a random position where is depart 
from the outer radius of the growing pattern. Then the 
particle jumps along the direction from the current site 
to one of its linked sites which are not occupied by the 
growing pattern, with equal probability step by step. At 
last, the particle irreversibly sticks at the nearest neigh- 
bor site of the growing pattern in terms of the physical 
distance and the pattern will grow gradually. To reduce 
the effect of fluctuation, the calculated result is taken 
average over 10 different network realizations and 10 in- 
dependent runs for each network configuration for each 
set of parameters {9,m). 



III. SIMULATION RESULT AND DISCUSSION 

Fig 2 illustrates the patterns of DLA which exhibit rich 
behaviors for different parameters 9 and m. For each 9, 
it can be seen that with the increase of m the patterns of 
DLA become thicker and denser, however, which is not 
obvious for large 9, approximately 1. For each m, the 
pattern is nearly the most dense when 9=0. While, it 
gets thin and sparse when 9 departs from tending to 1 
or —1. However, it is astonishing that the symmetry of 
the aggregation pattern is markedly broken while 9 < 0, 
which is more obvious as 9 tends to —1. To quantify 
the patterns of DLA, we calculate the fractal dimensions 
Dq of the DLA by box-counting method0,Q, which are 
shown in Fig. 3. It is clear that there exists a maximum 
value of Dq when 9 equals for each m. It can be seen 
that the more 9 departs from 0, the more Dq decreases. 
Moreover, it is found that Dq decreases more fast when 
9^1 than 9 -1. 

It is well-known that the special randomly branching, 
open structure of a DLA pattern results from the effects 
of screening, which is manifested through the fact that 
the tips of most advanced branches capture the incoming 
diffusing particles most effectively. In the present work, 
due to the long-range connections, particles can jump 
directly to distant sites, including the traditionally com- 
pletely screened deep fjord. The nonlocal connections ef- 
fectively weaken the screening effect so that the pattern 
of aggregate becomes compact and the fractal dimension 
Do increases with m increasing. On the other hand, the 
clustering exponent 9, which restricts the spatial length 
of the long-range connection, affects the morphology of 
the aggregate and the fractal dimension Dq. As 1, 
the long-range links are restricted more and more local in 
its vicinity so that the capacity of weakening the screen- 
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ing effect vanisfies graduafly. FinaUy, Dq does not vary 
and the morphofogy of the aggregate seems like the orig- 
inal DLA pattern, as shown in Fig. 2 and Fig. 3. When 
9=0, the spatial lengths of the long-range connections 
are entire random, neither too distant nor too local and 
they have a uniform distribution. Then, due to intensive 
weakening for the screening effect, the random particle 
has the chance to appear on arbitrary sites on the un- 
derlying network so that the pattern becomes thick and 
compact, corresponding to increase of Dq. However, as 
6^—1, the long-range links tend to the sites as distant 
as possible and the irreversible jumps along directed links 
break the symmetry of dynamics. Thus, small fluctua- 
tions are enhanced, and this instability together with the 
randomness inherent in the model leads to a complex 
asymmetrical behavior. 

However, the fractal dimension Dq is a rough de- 
scription because the pattern becomes asymmetric while 
9 —>■ —1. So, we have performed the multifractal 
analvsepM 0| here to see more details. It should be 
noted that our measurements concern the pattern itself 
other than its harmonic measure 25). 

Further characterization of the DLA patterns can be 
achieved by determining the generalized fractal dimen- 
sions Dq. Cover the pattern with a grid of square boxes 
of size e and define Pi{e) to be the relative portion of the 
pattern in cell i, and define N{e) to be the total number 
of boxes of size e needed to cover the whole pattern. The 
relative portion Pi{e) can be described as multifractal as: 



where a is the singularity, Na{s) the number of small 
squares of relative size e with the same singularity, and 
/(a) is the fractal dimension. 

To describe quantitatively the subtle geometrical fea- 
ture of the pattern, the fractal dimension /(a) can be 
obtain from the partition function Xq- 



N{e) 



Xq 



and its power law of £, 



In 



Xq 



Tq = Hm 

e~*o \ne 



(4) 



(5) 



where q is the moment order and Tq the index of the pow 
law. The generalized fractal dimension is defined as: 



D -— ^ 



(6) 



Then, (a, f{a)) can be obtained from (q, Dq) by perform- 
ing the Legendre transformation: 



a = -[{q - l)Dq], 



f{a) = aq- {q- l)Dq, 



(7) 



(8) 



In our calculation, the moment order q is taken for -30 

(2) to 30. 

We have calculated the multifractal spectra f{a) of 
the distribution of the patterns, with various number of 
long-range connections, m=1000, 10000, 100000, respec- 

(3) tively, for a original directed 2D SW lattice, 9=0. Figure 
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TABLE II: Some multifractal parameters of Figure 5. 



-1 -0.01 0.01 1 



^max 
A/ = f{amin) 



1.532 1.630 1.726 1.648 1.392 
2.969 3.307 3.591 3.109 2.727 
1.437 1.677 1.865 1.461 1.335 
0.352 0.652 0.752 0.676 0.347 
0.438 0.326 0.180 0.635 0.470 
fiamax) -0.086 0.326 0.572 0.041 -0.123 



IV. CONCLUSION 



a 



FIG. 5: The multifractal spectra with several value of the 
cluster exponent 9, -1, -0.01, 0, 0.01, 1, for m^lOOOO. 

4 shows the result. It can be seen that the curve becomes 
higher and the range of singularity a becomes wider with 
increasing the number of the long-range connections. In 
Table I, the multifractal parameters of the distribution 
are listed. The multifractal spectrum can be used to pro- 
vide more information about the subtle geometrical dif- 
ference, because of the amax and amin connecting with 
the smallest probability and the largest probability of the 
spatial distribution [show by Eq .(2)]. The result (Table 
I) show that amax and amin both increase with increas- 
ing m, while Aa=amax — amin also increases, indicating 
that the pattern becomes less irregular, less nonuniform, 
and more dense. Moreover, Fig. 5 illustrates that the 
multifractal curves with several value of the cluster ex- 
ponent e, -1, -0.01, 0, 0.01, 1, for a directed 2D SW lat- 
tice, TO=10000. More quantitative details can be seen in 
Table II. It illustrates that the range of a is the broadest 
and the curve is the maximal when 6=0, suggesting the 
pattern is the most dense, compact and regular, which 
corresponds to Fig. 2 and Fig. 3 showed above. 



TABLE I: Some multifractal parameters of Figure 4. 



m 




1000 


10000 100000 






1.598 


1.726 1.786 


^max 




3.332 


3.591 4.034 






1.734 


1.865 2.248 


f iS^min) 




0.684 


0.752 1.251 


f iS^max^ 




0.303 


0.180 0.147 


Af = f{a,nin) - 


f i^C^rnax , 


0.381 


0.572 1.104 
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